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Предложен метод краткосрочного прогнозирования с помощью вариационного сглаживающего сплайна, позволяющий опери- 
ровать большим числом настраиваемых параметров. Исследована зависимость погрешности прогнозирования от значения па- 
раметров и найдены их оптимальные значения. Проведен сравнительный анализ метода с аппаратом искусственных нейронных 
сетей в задаче краткосрочного прогнозирования. 


Постановка задачи 

Во многих задачах обработки измерений, таких, 
как интерпретация физических экспериментов или 
анализ статистических данных, возникает необхо- 
димость исследовать полученную зависимость с це- 
лью выявления характерных признаков и постро- 
ения прогнозирующей кривой. В настоящее время 
существует немало различных методов прогнозиро- 
вания, основные из них подробно рассмотрены в 
[1]. Среди математических методов прогнозирова- 
ния выделяют методы: экстраполяции, экспонен- 
циального сглаживания, прогнозирования по оги- 
бающей кривой, производственных функций, ис- 
следования операций, множественного регресси- 
онного анализа, сетевые методы, а также экономе- 
трические модели. Одним из отличий предложен- 
ного в [2] подхода к прогнозированию стало ис- 
пользование аппроксимирующих кусочно-полино- 
миальных функций, которые заменили регрессион- 
ные соотношения при моделировании и анализе 
процессов, что весьма важно для поиска экстрему- 
мов, прогнозирования и других задач обработки 
информации. Из всех кусочно-полиномиальных 
функций предпочтение отдается сплайн-функ- 
циям. Отличительная особенность сплайнов за- 
ключается в том, что они состоят из отрезков сте- 
пенного полинома малого порядка (степени), кото- 
рые сопряжены (склеены) в заданных узловых точ- 
ках процесса (узлах решётчатой функции). 

Для формального (трендового, детерминиро- 
ванного) прогнозирования важен тип исследуемого 
процесса. Например, при одних внешних условий 
процесс показывает экспоненциальный рост, а при 
смене их динамика показателя становится линей- 
но-падающей, что обнаруживает существенную не- 
стационарность. Чтобы собрать статистику поведе- 
ния показателей и продолжить её в прогнозном по- 
строении, требуется сделать отчётный период до- 
статочно долгим. Однако на длинном периоде уве- 
личивается вероятность изменений внешних усло- 
вий, он неизбежно делится на отдельные «куски», 
внутри которых поведение процесса может карди- 
нально меняться. Поэтому особую актуальность в 
моделировании, анализе и прогнозировании прио- 
бретает выбор единой (унифицированной, универ- 
сальной) системы приближающих функций. 


В настоящей работе рассмотрены существующие 
методы прогнозирования, проведен сравнительный 
анализ эффективности применения нейронных се- 
тей и вариационного сглаживающего сплайна 
(ВСС) при решении задачи прогнозирования. Также 
исследованы возможные режимы прогнозирования 
в реальном масштабе времени на основе ВСС для 
различных входных модельных данных и найдены 
оптимальные значения параметров ВСС. 

Расчетная схема и режимы функционирования 

вариационного сглаживающего сплайна 

Исследования, проведенные в [2], показали хо- 
рошую адаптивность сплайнов к моделированию, 
анализу и прогнозированию различных процессов. 
Унифицированный аппарат сплайн-функций за- 
менил широкий спектр аппроксимирующих мно- 
гочленов. Сплайны на интервалах имеют инвари- 
антную внутреннюю структуру, универсальную ма- 
тематическую форму, они самоподобны (как фрак- 
талы), поэтому удачно аналитически преобразуют- 
ся и реализуются, образуя после сопряжения еди- 
ный аппроксимирующий ансамбль. 

Предложенный в данной работе подход предпо- 
лагает построение ВСС на основе данных, посту- 
пающих в реальном масштабе времени и после- 
дующее вычисление прогнозируемых значений. В 
работе использовалась рекуррентная процедура 
построения ВСС для обработки текущих измере- 
ний, изложенная в [3]. Экстремальный функцио- 
нал А.Н. Тихонова /(5) для предлагаемого сплайна 
имеет вид [4]: 

/(У)=(і - р)(мо 2 |[у"(0] 2 л+рХ^(Ф-Ж)] 2 , 

,< /= о 

ч> 

где р - весовой множитель, устанавливающий 
компромисс между сглаживающими и интерполя- 
ционными свойствами сплайна 8(і), ре [0,1); А Г - 

р -р 

интервал дискретизации измерений; И = и 0 - 

число измерений внутри /-го звена сплайна, в даль- 
нейшем /г=соп$1, для любого /; Ц, Ц - начало и ко- 
нец /-го звена соответственно; у( //)=/(()+ <§(//) - из- 
мерения, представляющие собой смесь полезного 
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низкочастотного сигнала /(/) и широкополосной 
помехи ф). В задаче прогноза возможны исследо- 
вания и в условиях отсутствия помехи. Вариацион- 
ный подход и дал название сплайну. 

Значение сплайна в любой внутренней точке т 
для /-го звена определяется 

5Д т ) = а‘ 0 + а[т + а‘ 2 т 2 + а! ъ т ъ , 

О < т < к - 1 . 

Разрывные коэффициенты сплайна находятся 
из условия минимума целевого функционала в за- 
висимости от дефекта сплайна сі, то есть числа раз- 
рывных производных. Для дефекта й = 2 разрывны- 
ми являются коэффициенты а 2 ' и а' ь для дефекта 
<1 = 1 соответственно только а\. Непрерывные коэф- 
фициенты находятся из условия сопряжения 
звеньев сплайна. Расчетные формулы для всех ко- 
эффициентов приведены в [3]. 

Другим параметром ВСС является момент со- 
пряжения звеньев сплайна ц. Возможно несколько 
режимов использования ВСС в зависимости от со- 
отношения параметров дит. В данной работе ис- 
пользуется только текущий режим функциониро- 
вания ВСС, при котором вычисление ВСС осу- 
ществляется при поступлении каждого измерения 
с использованием (к-\) предыдущих измерений, 
сопряжение при этом осуществляется в момент 
времени і^\=Ц, где ц - номер дискретного отсчета 
і - го звена, в котором осуществляется сопряжение 
непрерывных производных сплайна. Схематично 
такой режим изображен на рис. 1. 



Рис. 1. Текущий режим ВСС: * - измерения, • - вычисленные 
значения, □ - момент д сопряжения звеньев ВСС 


Запаздывание результата Ѳ при обработке ин- 
формации для момента времени т {і‘ т ) складывает- 
ся из двух слагаемых: собственного времени вычи- 
сления оценок коэффициентов и значений ВСС Ѳ :> 
и времени запаздывания вычисления по отноше- 
нию к текущему времени, т. е. Ѳ=Ѳ 0 +(к-т)Аі. 

Для получения прогнозных значений использо- 
вание ВСС имеет ряд особенностей. Интуитивно 


ясно, что число измерений звена сплайна к - один 
из главных параметров настройки ВСС. Пусть 
необходимо спрогнозировать процесс вперед на і 
шагов квантования при известных к—і измерениях. 
Однако, алгоритм вычисления коэффициентов 
ВСС требует, чтобы все измерения в звене были из- 
вестны. Для решения возникшей задачи сократим 
интервал вычисления до к-і, используя главную 
особенность ВСС — возможность сопряжения 
звеньев сплайна в любой внутренней точке. При 
этом интервал построения сплайна остается преж- 
ним. Таким образом, при построении прогнозируе- 
мой кривой ВСС строится на всем интервале к- АТ, а 
коэффициенты вычисляются только для интервала 
{к~і)-Аі. Данный способ вычисления схематично 
изображен на рис. 2. 



Рис. 2. Построение кривой прогноза: А - прогнозируемые 
значения 


Влияние параметров вариационного сглаживающего 
сплайна на качество прогноза 

Используемый для прогнозирования ВСС име- 
ет несколько параметров настройки, существенно 
определяющих его свойства. Этими параметрами 
являются: сглаживающий множитель р, число из- 
мерений звена сплайна к, число прогнозируемых 
точек Определим некоторый набор оптимальных 
значений указанных параметров и полученное при 
этом качество прогноза. 

В работе для исследования параметров ВСС бы- 
ли выбраны две модельные функции: 

/, = 8 + 8, 4 ■ / - 0, 2 • г 2 + 1, 33 • 1 0 ' 3 • г 3 , (1) 


/, 


10-5ІП 



( 2 ) 


Для оценки качества прогнозирующих свойств 
сплайна использована приведенная среднеквадра- 
тическая погрешность, выраженная в процентах 
ст пр , %. Данная погрешность определяет меру бли- 
зости между модельными функциями (1) или (2) и 
результатом прогноза. 
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В работе проведен анализ влияния параметров 
ВСС на погрешность прогнозирования. Измене- 
ния значения ошибки прогнозирования при раз- 
личных р показаны на рис. 3 для функции (1). 



Прогноз при этом составляет три точки (г=3) 
при /г=10. Как видно, качество прогноза повыша- 
ется с увеличением р, т. е. при переходе ВСС в ин- 
терполяционный сплайн. Минимальное значение 
ошибки при этом достигает 0,48 %. Для функции 
(2) кривая зависимости имеет тот же вид, но значе- 
ния ошибки выше 1,19... 11,25 %. Таким образом, 
для периодической входной зависимости погреш- 
ность увеличивается примерно в два раза. 

Интересная закономерность наблюдается при 
изменении интервала к при постоянном числе то- 
чек прогнозирования. Оптимальное значение дан- 
ного параметра всегда составляет удвоенное коли- 
чество точек прогноза, т. е. непосредственно зави- 
сит от параметра і. Данное свойство можно наблю- 
дать на рис. 4 при г=3. 



Рис. 4. Влияние параметра к на ошибку 


Видно, что функция зависимости имеет один 
экстремум при Н= 2т, При ре[0,6;0,8] ошибка нес- 
колько снижается при увеличении /г, однако эк- 
стремум остается. Результаты аналогичны для об- 
еих функций, при этом погрешность для функции 
(2) также выше в два раза по сравнению с функци- 
ей (1). Для других соотношений к и жданная зако- 
номерность сохраняется. Возможно, это следует из 
теоремы Котельникова, так как анализируемые 
данные дискретизированы. Однако в других зада- 
чах анализа информации (фильтрация, сглажива- 
ние) с использованием ВВС подобных эффектов не 
наблюдалось [3]. 


Вообще говоря, выбор параметров ВСС осно- 
вывается на требуемом соотношении качество 
прогноза/длительность прогноза. Увеличивая чи- 
сло точек прогноза, мы тем самым снижаем его ка- 
чество. Однако, предложенный нами способ по- 
зволяет получить высокое качество прогноза при 
достаточно широком выборе количества прогнози- 
руемых значений. Пример зависимости погрешно- 
сти от длительности прогноза при /7=10 и р=0,9 для 
функции ( 1 ) показан на рис. 5. 



1 2 3 4 5 6 

Рис. 5. Влияние параметра г на ошибку 

Как видно, при прогнозировании первых трех 
точек погрешность держится примерно на одном 
уровне и составляет всего около 0,35 %. Для перио- 
дической функции - около 0,7 %. 

При проведении данного исследования также 
было выявлено, что при г, меньшем или равном 
половине интервала й, оптимальное значение сгла- 
живающего множителя составляют 0,6. ..0,8, т. е. 
уже оказывают существенное влияние сглаживаю- 
щие свойства ВСС. 



25 30 35 40 45 


Рис. 6 . Построение прогнозируемой кривой: - — - - вход- 
ная функция; в-н-н- прогнозируемые значения 

Добавим к задаче прогнозирования задачу сгла- 
живания входного сигнала (как известно, ВСС реша- 
ет ее весьма эффективно), в том числе и в реальном 
масштабе времени. Для этого внесем в исходные 
данные помеху, составляющей 10 % от максимально- 
го значения функции. С увеличением длительности 
прогноза погрешность прогнозирования увеличива- 
ется: при изменении гот 1 до 6 ошибка изменяется от 
1,8 % до 4 % при оптимально подобранных параме- 
трах. Таким образом, при прогнозировании трех то- 
чек ошибка составляет всего 2,2 %, что является не- 
плохим результатом, учитывая, что входной сигнал 
искажен помехой. Полученные значения указаны 
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для функции (1). Для функции (2) значение погреш- 
ности прогнозирования в среднем в два раза больше 
при аналогичных значениях параметров. 

Для наглядного представления построения 
прогнозируемой кривой приведем часть экспери- 
ментальных данных (рис. 6). 

На данном рисунке представлен фрагмент 
функции (1), не искаженной помехой, и прогнози- 
руемые значения на каждом интервале, й= 10, ^=3. 

Анализ эффективности нейронных сетей 

и ВСС при решении задачи прогнозирования 

Сопоставим эффективность ВСС и искусствен- 
ных нейронных сетей в задаче краткосрочного прог- 
нозирования. Построение нейронных сетей прово- 
дилось в программе-эмуляторе №ига1 №і\ѵогк \Ѵі- 
2аг<3 (ЫІЧ\У). В программе реализована многослой- 
ная нейронная сеть, обучаемая по алгоритму обрат- 
ного распространения ошибки. На базе N14^ были 
разработаны некоторые системы прогнозирования. 

При работе с нейронной сетью самым важным эта- 
пом является сбор статистики по процессу и обработ- 
ка входных данных. Обычно стараются собрать как 
можно больше данных, в какой-то мере влияющих на 
исследуемый процесс, а затем, в случае необходимо- 
сти, отсеять ненужные. Однако в случае краткосроч- 
ного прогноза набор входных данных невелик, чаще 
всего это несколько временных значений. 

Пусть необходимо спрогнозировать зависи- 
мость, подобную функции (1). На вход сети мы мо- 
жем подать только текущий момент времени, ха- 
рактеристики процесса и условия его протекания 
нам неизвестны. Создадим выборку из первых ста 
значений функции, в качестве обучающей выборки 
используем 50 % данных из нее, остальные данные 
будут тестовыми. На выходе необходимо получить 
значение функции в некоторый момент времени. 
Эксперимент показал, что сеть может найти значе- 
ние функции с точность до 0,5 % в 98 % случаев. 
Однако сеть работает только в границах обучаю- 
щей выборки, «в будущее» она смотреть не способ- 
на. Это является существенным ограничением 
нейронной сети, несмотря на то, что с задачей ин- 
терполяции она справилась хорошо, восстановив 
50 % промежуточных значений функции. 
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Если же провести подобный эксперимент над 
синусоидальной функцией (2), то в качестве об- 
учающей выборки необходимо задать более поло- 
вины имеющихся значений, иначе сеть не обнару- 
жит пики синусоиды. 

Вообще, как показал эксперимент, чтобы ней- 
ронная сеть прогнозировала будущее, ей необходи- 
мо большое количество информации, как в про- 
цессе обучения, так и в процессе работы, т. е. она 
применима для долгосрочного прогнозирования, 
когда имеется достаточная статистика по процессу. 
Еіапример, если мы узнаем период синусоиды и 
подадим его на вход сети, она сможет и в будущем 
распознать данную функцию. 

Таким образом, главное отличие ВСС и нейрон- 
ной сети заключается в представлении входных 
данных. ВСС требует только последовательные на- 
бор измеренных значений сигнала, и по ним, в си- 
лу гладкости входной функции, может построить 
участок функции «в будущем». При этом ВСС не 
сможет обнаружить периодичность функции, что 
также доказал наш эксперимент: в случае синус- 
оиды результаты прогнозирования оказались суще- 
ственно хуже. 

Нейронная сеть требует характеристики про- 
цесса, зависящие от периода, а не конкретного зна- 
чения времени. Она хорошо применима для перио- 
дических процессов, если входные данные содер- 
жат информацию о периодичности. Только в этом 
случае, нейронная сеть распознает синусоиду (2). 

В заключение отметим, что области эффектив- 
ного применения ВСС и искусственных нейрон- 
ных сетей существенно различны. Нейронные сети 
предъявляют высокие требования к входным дан- 
ным, необходимы большие затраты времени и ре- 
сурсов на их обучение. При этом нейронная сеть 
будет работать только в пределах обучающей вы- 
борки, т. е. будет применима лишь для того процес- 
са, на котором была обучена. Использование ВСС 
для прогнозирования универсально, поэтому дан- 
ный метод может быть применен в любой области, 
для любого процесса при минимальном наборе 
входных данных и без предварительных затрат вре- 
мени. Однако в силу своей природы, ВСС более 
эффективен для прогнозирования нестационарных 
процессов, чем периодических функций. 
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